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Nonparallel Stability of Boundary Layers 

The asymptotic formulations of the nonparallel linear stability of 
incompressible growing boundary layers are critically reviewed. These 
formulations can be divided into two approaches. The first approach 
combines a numerical method with either the method of multiple scales, 
or the method of averaging, or the Wentzel-Kramers-Brillouin (WKB) 
approximation; all these methods yield the same result. The second 
approach combines a multi-structure theory with the method of multiple 
scales. Proponents of the second approach have claimed that their 
approach is rational and the first approach is not rational. The first 
approach yields results that are in excellent agreement with all 
available experimental data, including the growth rates as well as the 
neutral stability curve. On the other hand, the second approach cannot 
even yield the neutral curve for the Blasius flow. 

Introduction 


This paper addresses the linear stability of incompressible growing 
boundary layers. For two-dimensional mean flows, the streamwise 
velocity component U(x,y) is a function of the transverse coordinate y 
as well as the streamwise coordinate x. However, the rate of variation 
of U with respect to x (i.e., all/ax) is small compared with the rate of 
variation of U with respect to y (i.e., all/ay). Moreover, the transverse 
velocity component V(x,y) is small compared with U and is a function of 
y as well as x. For three-dimensional flows, the velocity components 
U(x,y,z) and W(x,y,z) in the plane of the body are much larger than the 
transverse velocity component V(x,y,z). Moreover, all/ax, aU/az, aW/ax, 
and aW/az are small compared with aU/ay and aW/ay. 


To determine the linear stability of a three-dimensional mean flow, 
we superimpose on it a small disturbance u(x,y,z,t) v(x,y,z,t), 
w(x,y,z,t), and p(x,y,z,t). Substituting the total flow into the 
Navier-Stokes equations, subtracting the mean-flow quantities, and 
linearizing the resulting equations, we obtain 
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where velocities, lengths, and time were made dimensionless using the 
free-stream velocity U , a characteristic length a, and a characteristic 
time 6/U . Here, R = iTfi/v is the Reynolds number. The boundary 
conditions are 


u = v= w = 0aty = 0 (5) 


u, v, w, p -*■ 0 as y -*• ® (6) 

The terms in the square brackets in Eqs. (2)-(4) are due to the growth 
of the boundary layer (nonparallel terms). 

Parallel Problem 


Considering the parallel problem, one neglects the terms in square 
brackets in Eqs. (2)-(4) and considers U and W to be functions of y 
only. Then, one seeks a normal mode solution of the form 


u = ^(yJE, v = c 3 (y)E, w = e 5 (y)E, p = t 4 (y)E (7) 

where 

E = exp[i(ax + bz - wt) ] (8) 

Substituting Eqs. (7) and (8) into Eqs. (l)-(6) and neglecting the terms 
in square brackets yields 

iotCj + 0c 3 + iBs 5 = 0 (9) 

i (all + BW - u)c l + C 3 DU + ia? 4 - ^ (D 2 - a - B 2 )c t = 0 (10) 

i(aU + BW - u>)c 3 + Dc, - i (D 2 - a' - fi')^ = 0 (11) 

i (all + BW - w)c 5 + 5 3 DW + i6c 4 - ^ (D 2 - a - 8 )c 5 = 0 (12) 

Cx = C 3 = c 5 = 0 at y = 0 (13) 

c -*■ 0 as y -*• ® (14) 


where D = 3 /sy. For a given U(y) and W(y), Eqs. (9) - (14) constitute an 
eigenvalue problem, which yields a dispersion relation of the form 
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id = u>(a,6,R) 


(15) 

A number of techniques have been developed for solving this eigenvalue 
problem. These include shooting techniques, finite-difference methods, 
Galerkin methods, and collocation techniques using Chebyshev or Jacobi 
polynomials. 

For the case of a two-dimensional mean flow and a two-dimensional 
disturbance, W = 0, c 5 = 0 and 8 = 0, and Eqs. (9)- ( 14) reduce to 


iotCj + Dc 3 = 0 (16) 

i(aU - a»)c l + CjDU + iac 4 - I (D 2 - a 2 ) C 3 = 0 (17) 
i(aU - iu)c 3 + Dc 4 - | (D 2 - a 2 )c 3 = 0 (18) 
Cj = c 3 = 0 at y = 0 (19) 
G n - 0 as y - » (20) 


Equations (16)-(20) can be combined to yield the Orr-Sommerfeld equation 
(i a R) _1 (D 2 - oi 2 ) 2 ^ 3 = (U - c)(D 2 - a 2 )c 3 - c 3 D 2 U 


( 21 ) 

subject to the boundary conditions 

C 3 = 0c 3 = 0 at y = 0 (22) 


c 3 , Dc 3 - 0 as y » (23) 

where c = u>/a. The neutral stability curve calculated using either Eqs. 
(16) - (20) or Eqs. (21)-(23) are in good agreement with available 
experimental data as shown in Figure 1. 

Recently, Smith^ claimed the above methods to be "irrational" and 
developed multi-structured theories for treating this problem. He used 
a result from an "irrational theory" for the Blasius flow to observe 
that "the typical wavelength of the neutrally stable modes on the lower 
branch increases proportionally to Re ' as Re - and concluded that 
disturbances at the 10^ branch vary on a streamwise length 0(Re -3 /®) 
and a time scale 0(Re" ' ) and hence they are governed by a triple-deck 
structure. Consequently, he let 

x = 1 + e 3 X, t = e 2 T (24) 


and 
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u,v,p<*E = exp[ie(x) - inT] 

(25) 

where 

e = Re -1 ^ 8 and 



2 3 



ft = oj 1 + eu> 2 + e u> 3 + e line a> 4 ^ + 

(26) 


de 2 3 



^ = a x + ea 2 + e a 3 + e fi,ne a 4 ^ + 

(27) 

Then, 

he expanded the variables in the three decks as follows: 


Main Deck 



u = [uj + eii 2 + e 2 u 3 + e 3 jine u 4L + ...]E 

(28) 


r 2 3 4 i.- 

V = [eV 1 + e V 2 + e V 3 + e fcne V 4 ^ + ...]E 

(29) 


p = [ep 1 + e 2 p 2 + e 3 p 3 + e\ne p 4( _ + ...]E 

(30) 

where 




y = e 4 Y, Y = 0(1) 


Lower 

Deck 



U = [ll 1 + e: U 2 + e U 3 + e Jine U 4 ^ + *..]E 

(31) 


V = [ e V : + c M 2 + e V 3 + e line V 4 ^ + ...]E 

(32) 


P = l e ^i + e P 2 + e P 3 + e &ne P 4 ^ + ...]E 

(33) 

where 




y = e 5 Z, Z = 0(1) 


Upper 

Deck 



r - 2_ 3_ 4 , 

U - [eU 1 + e U 2 + e U 3 + e £ne U 4 ^ + .•.]£ 

(34) 


r - 2_ 3_ 4 , 

V - [eV l + e V 2 -f e V 3 + e line V 4 ^ + ...JE 

(35) 


p = [ep x + e p 2 + e p 3 + e ine P 4 ^ + ...]E 

(36) 


where 


y = e 3 y. y = o(i) 
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To account for the nonparallel effects. Smith had to include the next 
term in each of the expansions in Eqs. (26)-(36). 

Substituting the above expansions into the parallel part of the 
disturbance equations ( 1) -(4) and boundary conditions (5) and (6), 
dropping the terms in square brackets, putting W = 0, separating 
coefficients of like powers of e, and solving the resulting 36 
equations, one obtains expressions for u, v, and p in the different 
decks. Matching the resulting expressions provides asymptotic 
expansions for u, v, and p. For the neutral stability curve. Smith 
obtained 

F = lu = 0.995 R ~y 2 [1 + 1.597 R“ ^ + 10.02 R“ ** 
n n 6 1 6 6 

+ 0.988 R in R + ...1 
6 6 1 


where 


R = 1.7208/xRe 

O 

This expression is in fair agreement with the lower branch of the 
neutral stability curve for large Re. However, its accuracy 
deteriorates as Re decreases. In fact, it does not predict a minimum 
critical Reynolds number. 

Bodonyi and Smith^ inspected the results of the "irrational theory" 
to observe that "the stability properties of the Blasius boundary layer 
are ggv^rned by the behavior on the streamwise length scale 
0(Re" ' ) as far as the upper branch of the neutral curve is 

concerned". Consequently, they used o = Re~' as their perturbation 
parameter and us|d the streamwise scale X defined by x = 1 + o 9 X and the 
time scale t = o T. This choice leads to a five-zoned structure. To 
account for the nonparallel effects, one needs to carry out the 
expansion to 0(o ). In view of the logarithmic terms, one needs 13 terms 
in the expansion. With three variables and five decks, one needs to 
derive and solve 195 equations and then match the results. Bodonyi and 
Smith gave up after four terms. Their calculated neutral stability 
curve, which is intended to approximate the upper branch, is below the 
lower branch!! We note that for an Re = 10 6 , o * 0.5, which is not 
smal 1 . 

For the case of an accelerating b^ndary layer. Smith and Bodonyi^ 
assumed a streamwise variation 0(Re -5 ' ) and a time scale 0(Re“ ' ) 

near the upper branch of the neutral stability curve. Using this 
streamwise variation leads to a five-zoned structure, with the 
nonparallel effects appearing at 0(Re -5 ' ). 

It should be noted that the parallel flow assumption breaks down 
miserably for the case of Gortler instability. Floryan and Saric^ and 
Ragab and Nayfeh 5 derived the appropriate equations for Gortler 
instability for the cases of zero and nonzero pressure gradients, 
respectively. Hall® questioned the solution of the resulting equations 


2SO 



using a normal-mode approach and suggested solving them as an initial- 
value problem. 


Nonparallel Problem 

A better agreement between the theoretical and experimental results 
can be obtained by accounting for the influence of the nonparallel 
terms' -10 . To this end, we can use either the method of averaging or 
the WKB approximation or the method of multiple scales 11,12 . In this 
paper, we use the method of averaging and let 

u = A(x,z,t)^(y,x)e 19 , v = A(x,z,tH 3 (y,x)e 10 (37) 


p = A(x,z,t)c 4 (y,x)e 19 , w = A(x,z,tH 5 (y,x)e 19 
where A is a slowly varying function of x and t, 

w " a(x) ’ « “ 8 ’ It " ' “ 


(38) 


(39) 


and the c are given by Eqs. I9)-(14). After a straightforward but 
iTgebra, one obtains 1 ^ -1 ' 
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where the x, are solutions of the adjoint homogeneous problem. 
(40) can be n rewritten as 13 


(44) 

Equation 
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( 45 ) 
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where 



Here, and w. are the components of the group velocity in the 
streamwise direction. 


Equation (45) describes the propagation of a wavepacket centered at 
the frequency uy and the wavenumbers a r and $ r , where the subscript r 
stands for the real part. Nayfeh 14 showed that for a physical problem, 
w and oj in Eq. (45) must be real. 

a 6 'i \ / 


For a monochromatic wave, aA/at = 0 and Eq. (45) reduces to 


Lii 
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3X 


3A - hfl 
a) — = nA 
8 3Z 


(46) 


For a physical problem, Nayfeh 14 showed that a> B /u> a must be real. For 
the case of a parallel mean flow, this condition reduces to da/dB being 
real, which was obtained by Nayfeh 14 and Cebici and Stewartson 1 ® using 
the saddle-point method. 


Two-Dimensional Mean Flows 

For the case of a monochromatic wave, aA/at = 0 and Eq. (40) yields 

A = A exp[ J (h 4 /h 2 )dx] (47) 

x 

0 

where A is a constant. Hence, 

h 4 

u = A Q c 1 (y,x)exp[i(Jadx - cot) + J (|p)dx] (48) 

Consequently, the growth rate 

a = Real Ifx ^ ln U ^1 

is given by 

h. 

« = Real (^) + Real [|j (ln t ,)| (49) 

The first term is the quasiparallel contribution, whereas the last two 
terms are due to nonparallelism. It should be noted that the last term 
produces a variaton in the growth rate across the boundary layer. 

Since Ci is a function of y and, in general, distorts with 
streamwise distance, one may term stable disturbances unstable or vice 
versa. Morever, a different growth rate would be obtained if one 
replaces u with another variable. For example, using v or p or w, one 
obtains the growth rates 

h 4 

° - - «, + R^' (hj) + Real [fj- (lncj] (50) 
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where m = 3, 4, and 5, respectively. This raises the questions "What is 
meant by stability of a boundary layer?" If the stability criterion is 
based on o, then which a should be used? If one uses an N factor to 
compare the stabilizing or destabilizing influences of certain 
modifications to the boundary layer, then the contribution of the last 
term will not be significant. 

In the case of parallel flows, the last terms in Eqs. (49) and (50) 
vanish and the growth rate is unique and independent of the variable 
being used. Consequently, one can speak of neutral disturbances or 
neutral stability curve given by the locus of a.j(R,w) = 0. However, in 
the case of a nonparallel flow, the neutral stability is given 
by a(R,w) = 0 and depends on the flow variable used to calculate the 
growth rate and the distance from the wall. To compare the analytical 
results with experimental data, one needs to make the calculations in 
the same manner in which the measurements are taken. Available 
experimental stability studies almost exclusively use hot-wire 
aneomometers. Usually, they measure the rms value |u| of the streamwise 
velocity component u and use it to define the growth rate. Figure 2 
compares the neutral stability curves calculated using |u| and + Q 2 + v 2 
with the experimental data of Kachanov, Kozlov and Levchenko 19 . .Since 
the experiment measured |u|, the calculations of Saric and Nayfeh 10 , 
which were based on |u[, are in better agreement with the experimental 
data than the calculations of Bouthier , which were based on Q 2 + v. 
Moreover, the growth rate is singular at the locations where |u| =0. 
Figure 2 shows also that the calculated locations of the singular growth 
rates are in good agreement with the experimental results. 

Some of the available experimental studies follow the maxima of |u| 
whereas others follow a constant bound ary- layer similarity variable n. 
Saric and Nayfeh 9 ’ 10 found that the contribution of the last terms in 
Eqs. (49) and (50) are significant if one follows a constant n whereas 
their contributions are negligible if one follows the maxima of |u|, 
yielding 


a = - a. + Real(^) (51) 

The neutral stability curve calculated by Saric and Nayfeh 9 using Eq. 
(51), and shown in Figure 1, is in very good agreement with the 
experimental data that follow the maxima of |uj, except near the minimum 
critical Reynolds number where the data may be suspect. However, in the 
case of experiments conducted by following trajectories of constant n 
such as those of Ross et al. , the effect of the distortion of the 
eigenfunction cancels the nonparallel effects, resulting in a better 
agreement between their data and the results of quasiparallel theory. 

Saric and Nayfeh 10 made other comparisons of the growth rates 
calculated using Eq. (51) with the experimental data of Strazisar, Prahl 
and Reshotko^ 1 and Kachanov, Kozlov and Levchenko 19 . Strazisar et al. 


* Also, private communication, June 1976. 
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conducted their experiments in a water tunnel and performed their 
measurements at the maxima of |u|, thereby minimizing the effects of the 
distortion of the eigenfunction. They measured the amplification rate 
as a function of the frequency at different locations on the plate, 
corresponding to different Reynolds numbers. Figure 3 shows a good 
agreement between the theoretical and experimental results. Kachanov et 
al . also followed the maxima of |u| and measured the amplification 
factor a = |u|/|uJ, where |u Q | is the rms value of u at the first 
neutral point. Figure 4 shows a good agreement between the theoretical 
results calculated by Saric and Nayfeh^" using Eq. (51) and the 
experimental results. 

The present nonparallel analysis was extended by El-Hady and 
Nayfeh^ to the case of two-dimensional compressible boundary layers, by 
Nayfeh 14 to the case of three-dimensional compressible boundary layers, 
and by Nayfeh and El-Hady^ and Asrar and Nayfeh^ to heated boundary 
layers. 
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Figure 1. Neutral stability curve for Blaslus boundary layer. Solid 
symbols are Branch I experimental points. Open symbols are 
Branch II. The critical Reynolds number is 400 for 
nonparallel calculations, 520 for parallel calculations 
(Saric and Nayfeh^). 
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Figure 2. Vertical variation of neutral stability points at F = 

200x10”®. Experimental points from Kachanov, Kozlov and 

Levchenko*^. Dashed lines are calculations of Bouthier^ 

based on energy. Solid lines are calculations of Saric and 

Nayfeh*® based on |u|. Streamwise position is the Reynolds 

number based on 5„ which is the s of Refs. 7 and 17. Solid 
r 

triangles give the locus of |u| = 0 and the broken line is 
the calculation*® for |u| = 0. 
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Figure 3. Amplification rate as a function of frequency: Theory of 

Saric and Nayfeh; Experimental data of Strazisar, Prahl 
and Reshotko following maximum of |uj. 
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Figure 4. Amplification factor a as a function of streamwise position 
at F = llOxlO" 6 and F = 86 x 1CT 6 . Experimental points from 
Kachanov, Kozlov and Levchenko following maximum of |u| and 
nonparallel results of Saric and Nayfeh based on following 
maximum of |u | . 
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